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The AMU Team completed one task, began another, and 
continued work on three: 


Launch Support 

Mr. Wheeler and Dr. Merceret 
supported the launch of Shuttle 
Discovery on 24 February. 

Dr. Bauman and Dr. Merceret 
supported the Atlas V launch 
on 5 March. 

Ms. Crawford and Dr. Merceret 
supported the Delta IV launch 
on 1 1 March. 


■ 


Dr. Watson completed the task to update the AMU-developed severe weather forecast tool, and be- 
gan work on the second phase of verifying the performance of the MesoNAM weather model over 
Kennedy Space Center and Cape Canaveral Air Force Station. 

Ms. Crawford continued work to improve the AMU peak wind tool by processing and analyzing wind 
tower data to determine peak wind behavior during times of onshore and offshore flow. 

Dr. Bauman continued updating lightning climatologies for airfields around central Florida and creat- 
ed new climatologies for specific moisture thresholds as defined by Florida soundings. 

Mr. Wheeler completed a study for the 30th Weather Squadron at Vandenberg Air Force Base in 
California in which he found precursors in weather observations that will help the forecasters deter- 
mine when they will get strong wind gusts at their northern towers. 
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Quarterly Task Summaries 

This section contains summarizes of the AMU activities for the second quarter of Fiscal Year 201 1 (January - 
March 2011). The accomplishments on each task are described in more detail in the body of the report starting 
on the page number next to the task name. 


Peak Wind Tool for User LCC, Phase IV 
( Page 5 ) 

Purpose : Recalculate the Phase III cool season peak wind statistics 
using stability as an added stratification. Peak winds are an important 
forecast element for launch vehicles, but the 45th Weather Squadron 
(45 WS) and Spaceflight Meteorology Group (SMG) indicate that 
they are challenging to forecast. Stability has long been known to 
have a strong affect on surface winds. Recalculating the statistics 
after stratifying by stability will make them more robust and useful to 
operations. 

Accomplished: Determined that the tower data cannot be used for 
the stability stratification. The Cape Canaveral Air Force Station 
(CCAFS) soundings are being processed to determine if they can be 
used to determine stability with more accuracy. The tower data were 
stratified by on/offshore flow to calculate hourly climatological values 
for the 5-minute mean and peak speeds, gust factors, and the num- 
ber of occurrences for each sensor on the towers and for each 
month. There appears to be a relation between the hourly values of a 
solar parameter and the gust factors. 




Situational Lightning Climatologies for 
Central Florida, Phase V ( Page 7 ) 

Purpose: Update the existing lightning climatology to improve op- 
erational weather support to Kennedy Space Center (KSC), 
CCAFS, Patrick Air Force Base (PAFB), and commercial and gen- 
eral aviation across central Florida. The update includes adding 
more years of data to the database, adding more sites and adding 
stratifications for moisture and stability parameters. These updates 
will provide climatologies for new sites for which the 45 WS and 
National Weather Service (NWS) have forecast responsibility, and 
to help forecasters distinguish lightning days that are more active 
from those that are less active within the same flow regime. 

Accomplished: Received National Lightning Detection Network 
(NLDN) data (May-September 1 989-201 0) from the 1 4 WS for the 
nine additional sites requested by the NWS in Melbourne, Fla. The 
data were stratified by precipitable water (PWAT) thresholds and 
new PWAT-based lightning climatologies were generated for 28 of 
the 32 airfields. 
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Quarterly Task Summaries 

(continued) 



Vandenberg Air Force Base North Base 

Wind Study ( Page 8 ) 

Purpose: Analyze local wind tower, surface, upper air and 
sounding data from Vandenberg Air Force Base (VAFB) to find 
precursors to high wind events in the north base towers. The 
30 WS states that terrain influences the unpredicted strong 
northeast winds that have been measured on several of the 
north base wind towers and exceed their 35 kt warning criteria. 
This study will examine those influences and document any 
precursors that may be found that will assist forecasters in an- 
alyzing their wind warning criteria. 

Accomplished: The VAFB wind tower data provided by the 
30 WS were decoded and evaluated. Thirty of the 66 event 
days had all the data needed for further analysis. Surface and 
upper air charts were analyzed to determine the synoptic con- 
ditions for each event day, and peak wind direction and speed 
charts were developed and analyzed. The results showed that 
all of the high wind events began with a wind shift to the north- 
east followed by the strong winds within 2 hours of the shift. 


Upgrade Summer Severe Weather Tool 
Phase III ( Page 11 ) 

Purpose: Upgrade the Summer Severe Weather Tool by add- 
ing another warm season and testing another statistical tech- 
nique to determine if its performance can be improved. This 
task increases the period of record from 21 to 22 years and 
uses logistic regression to determine the appropriate predic- 
tors and provide a probability forecast. The performance of the 
logistic regression equations will be compared with the previ- 
ous tool. 

Accomplished: Developed a four-predictor logistic regression 
equation by determining how each individual predictor contrib- 
uted to the reduction of variance. The equation’s performance 
was tested against that of the current severe weather forecast 
tool using four statistical tests. The results showed that the 
current tool outperforms the equation. 
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Quarterly Task Summaries 

(continued) 
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MesoNAM Verification Phase II 
( Page 13 ) 

Purpose : Update the current tool that provides objective verifi- 
cation statistics of the 12-km North American Mesoscale (NAM) 
model (MesoNAM) for CCAFS and KSC. This tool helps the 
Launch Weather Officers understand the model’s performance 
when they use it to evaluate launch commit criteria (LCC) dur- 
ing launch operations. The modifications include adding a year 
of observations and model output data to the original database. 
The objective analysis consists of comparing the MesoNAM 
forecast winds, temperature and moisture to the observed val- 
ues at the KSC/CCAFS wind towers used to evaluate LCC. 

Accomplished: The wind tower observations and MesoNAM 
forecasts needed for the task were acquired. The tower data 
were quality controlled (QC-d) and processed to remove un- 
needed time periods and fill in missing values. When the data 
were prepared, the mean value for each observed parameter 
was computed for each hour. 
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AMU ACCOMPLISHMENTS DURING THE PAST QUARTER 

The progress being made in each task is provided in this section, organized by topic, 
with the primary AMU point of contact given at the end of the task discussion. 

SHORT-TERM FORECAST IMPROVEMENT 


Peak Wind Tool for Us- 
er LCC, Phase IV 
(Ms. Crawford) 

The peak winds are an important 
forecast element for the Expendable 
Launch Vehicle and Space Shuttle 
programs. As defined in the Launch 
Commit Criteria (LCC) and Shuttle 
Weather Flight Rules, each vehicle 
has peak wind thresholds that cannot 
be exceeded in order to ensure safe 
launch and landing operations. The 
45th Weather Squadron (45 WS) and 
the Spaceflight Meteorology Group 
(SMG) indicate that peak winds are a 
challenging parameter to forecast, 
particularly in the cool season. To 
alleviate some of the difficulty in 
making this forecast, the AMU calcu- 
lated cool season wind climatologies 
and peak speed probabilities for 
each of the towers used to evaluate 
LCC (Figure 1) in Phase I (Lambert 
2002). In Phase III (Crawford 2010), 
the AMU updated these statistics 
with six more years of data, added a 
new time-period stratifications and 
created a graphical user interface 
(GUI) to display the desired values 
similar to that developed for SMG in 
Phase II (Lambert 2003). Based on 
recommendations from Phase III and 
observations by the launch weather 
officers (LWOs), the 45 WS tasked 
the AMU to stratify the data by stabil- 
ity and onshore/offshore flow and 
recalculate the climatologies and 
probabilities. These modifications will 
likely make the statistics more robust 
and useful to operations. 

Data 

Ms. Crawford received all Kenne- 
dy Space Center (KSC) and Cape 
Canaveral Air Force Station 
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(CCAFS) wind tower network and 
CCAFS sounding data for December 
2010 from Mr. Madison of Computer 
Sciences Raytheon (CSR) and the 
hourly Shuttle Landing Facility (SLF) 
pressure data from the 14 WS need- 
ed for the Richardson number calcu- 
lations. She quality controlled and 
prepared these data for analysis us- 
ing the S-PLUS (Insightful Corpora- 
tion 2007) statistical software pack- 
age. The data are from the cool sea- 
son months October through April 
1995-2010, 16 cool seasons in the 
period of record (POR). 

Stability Determination 

Ms. Crawford calculated the gra- 
dient and bulk Richardson numbers 
(Ri and R B ) for each level on Towers 
2, 6, 110, and 313 collected in Janu- 
ary, all years in the POR. Stull (1988) 
states that flow 
becomes turbu- 
lent when R, < 

0.25. Using this 
criterion, the re- 
sults showed the 
tower layers 
were unstable 
over 90% of the 
time at Towers 
2, 6 and 100, 
and 75% of the 
time at Tower 
313. The high 
percentages of 
unstable cases 
did not appear 
representative of 
the climatologi- 
cal stability for 
January. 

To examine 
these results fur- 
ther, Ms. Craw- 


ford determined the percentage of 
unstable values for each hour of the 
day at Towers 2 and 313, shown in 
Table 1 . There was a clear diurnal 
signal in the values. In January, the 
sun rises locally at 0700-0715 EST 
and sets at 1745-1800 EST. In 
Table 1, the night hours between 
sunset and sunrise are in the left 
three columns, and the day hours 
between sunrise and sunset are in 
the right three columns. The over- 
night values remained steady with a 
slight increase in the midnight hours 
between 1100 and 0300 EST in both 
towers, and values were 10-15% 
lower at Tower 313. The percent- 
ages increase quickly after sunrise to 
99-100% between 1000 and 1600 
EST, and then dropped quickly to 
overnight values by sunset. 


Table 1 . Hourly frequencies in percent of R, < 0.25 for 
Towers 2 and 313 in January 1995-2010. 


Hour 

EST (UTC) 

19 ( 00 ) 

20 ( 01 ) 

21 ( 02 ) 

22 ( 03 ) 

23 ( 04 ) 

00 ( 05 ) 

01 ( 06 ) 

02 ( 07 ) 

03 ( 08 ) 

04 ( 09 ) 

05 ( 10 ) 

06 ( 11 ) 


71 

70 

71 
71 

73 
76 
78 

76 

77 

74 


74 


74 


313 


57 

57 
59 
61 

63 
62 

64 
62 
62 

58 

57 

58 


Hour 

EST (UTC) 

07 ( 12 ) 

08 ( 13 ) 

09 ( 14 ) 

10 ( 15 ) 

11 ( 16 ) 
12 ( 17 ) 
13 ( 18 ) 
14 ( 19 ) 
15 ( 20 ) 
16 ( 21 ) 

17 ( 22 ) 

18 ( 23 ) 


73 

82 

97 

100 

100 

100 

100 

99 

99 

99 

94 

78 


313 


54 

63 

90 

99 

100 
100 
100 
100 
100 
100 
96 


64 






Figure 1. Surface weather map for 20 January 
2003 at 0700 EST (1200 UTC) from 

http: //www. hoc, ncep. noaa, gov/daiivwxmaD/ . 


The quick increase to 100% 
was likely due the surface-level 
instability created by solar heat- 
ing of the ground . The cause of 
the slight increase in the steady 
values around midnight is un- 
clear. Nonetheless, the percent- 
ages of unstable cases were still 
higher than expected for the 
overnight. Unstable R, values 
occurred both day and night 
even when the region was under 
the influence of a strong, stable, 
high pressure center (Figure 1). 
Based on these results, Ms. 

Crawford concluded that the tower 
data could not be used to determine 
stability and, therefore, could not be 
used to stratify the data by stability. 

It is possible that the CCAFS 
soundings can be used to determine 
boundary layer stability. Ms. Craw- 
ford sent the quality controlled (QC- 
d) soundings to Mr. Kienzle of EN- 
SCO, Inc.’s GeoSystem Solutions 
Division to calculate the mixed layer 
(ML; Stull 1988) height using algo- 
rithms he developed for transport 
and diffusion models. The ML height 
will serve as the proxy for the height 
of the boundary layer. 

On/Offshore Stratifications 

Ms. Crawford determined the 
times of onshore and offshore flow 
using each tower wind sensor 
(Figure 2), and then calculated the 
hourly means of the winds, gust fac- 
tors, and frequency for each flow pat- 
tern/month combination. Towers 2, 6, 
and 110 have two sensors at each 
height on opposing sides, NW and 
SE. The towers at pads 39A and B 
and Tower 41 have one sensor on 
each tower, but each pad has towers 
located on the NW and SE sides. 

This makes their configuration similar 
to Towers 2, 6, and 110. Ms. Craw- 
ford used the algorithms developed 
in Bauman (2010), which accounted 
for the sensor location relative to the 
tower and the wind directions: 

Tower 2; 

• Onshore if SE sensor direction is 

> 46° and < 225°; NW sensor di- 
rection used if SE sensor direc- 
tion missing. 


• Offshore if NW sensor direction is 

> 226° or < 45°; SE sensor direc- 
tion used if NW sensor direction 
missing. 

Towers 6, 1 10, 41, and 39X: 

• Onshore if NW sensor direction is 
£ 316° or < 45°, or SE sensor 
direction is > 46° and < 135°. 

• Offshore if NW sensor direction is 

> 226° and < 315°, or SE sensor 
direction is > 136° and < 225°. 

• If direction from one side missing, 
direction from other side used. 

Tower 108: 

• Onshore if direction is > 316° or 

< 135°. 

• Offshore if direction is > 1 36° or 

< 315°. 

The onshore/offshore flow criteria 
for Tower 2 are different than for the 
other towers due to its location near 
the south end of CCAFS where the 
coastline is oriented NE-SW, approx- 
imately 45° to 225° (Figure 2). The 
coastline nearest Towers 6, 108, 

110, 41 , and 39A and B are oriented 
NW-SE. approximately 315° to 135°. 
To increase the sample size, Ms. 
Crawford added the condition to use 
the direction from the downwind sen- 
sor if the direction from the upwind 
sensor is missing. There were some 
conflicts for towers with two sides in 
which one sensor indicated one flow 
type and the other sensor showed 
the opposite. This usually occurred 
during transitions between onshore 
and offshore flow. In these cases, 

Ms. Crawford assigned the most re- 
cent unambiguous flow pattern. 


Hourly Climatologies 

Figure 3 shows the hourly cli- 
matologies of the 5-minute av- 
erage and peak wind speeds 
(left) and the number of occur- 
rences of each flow pattern for 
onshore and offshore flow at 
the 90-ft NW sensor on Tower 
2 in January. The offshore av- 
erage and peak wind speeds 
show the expected increase in 
values during the daytime hours 
(1200-2300 UTC), but the on- 
shore values decreased by 1-3 
kt from their nighttime values 
during the same time period. The 
main feature in number of occur- 
rence curves is the marked increase 
in offshore flow events and decrease 
in onshore flow events beginning at 
midnight local time (0500 UTC) and 
peaking between 0400-0800 EST 
(0900-1300 UTC). This likely reflects 
the occurrence of the land breeze in 
the overnight hours known to occur 
over KSC/CCAFS. 

The gust factor means and 
standard deviations shown in Figure 
4 followed similar trends as the 
speed climatologies in Figure 3: rela- 
tively constant values during the 
night, increasing values with the sun- 



locations of the launch pads and 
LCC wind towers. 
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Hourly Mean 5-Min Average and Peak Winds 
On and Offshore Flow at Tower 2 NW 90 ft in January 
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Figure 3. Hourly climatological values for the wind speeds at the 90-ft NW sensor on Tower 2 in January: (left) the 
hourly mean 5-minute average and peak winds and (right) the hourly number of occurrence of onshore and offshore 
flow. Offshore values are blue and red, and onshore values are green and purple in each chart. 


rise to mid-day, and decreasing through sunset. The on- 
shore mean values were consistently less than the off- 
shore values, but the standard deviations were similar. 
Dr. Merceret compared the gust factor curves to a solar 
parameter for January. This solar parameter varies from 
0 to 1 and depends on the sun angle for each hour and 
day of year (http://www.qcstudio.com/suncalc.html ). 
Figure 4 also shows the solar parameter hourly curve 
calculated for 15 January (mid-month). There is a visual 
correlation between the GF and solar curves. Dr. Mer- 
ceret will continue to explore the relationship between 
the solar parameter and the gust factor means and 
standard deviations. 

Contact Ms. Crawford at 321-853-8130 or 

crawford.winnie@ensco.com for more information. 


Hourly Solar Parameter and Gust Factor 
Mean/Standard Deviation 

On and Offshore Flow at Tower 2 NW 90 ft in January 
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Figure 4. Hourly gust factor means and standard 
deviations at the 90-ft NW sensor on Tower 2 in 
January , and the solar parameter curve for January 15. 


Situational Lightning 
Climatologies for Cen- 
tral Florida, Phase V 
(Dr. Bauman) 

The threat of lightning is a daily 
concern during the warm season in 
Florida. Research has revealed dis- 
tinct spatial and temporal distribu- 
tions of lightning occurrence that are 
strongly influenced by large-scale 
atmospheric flow regimes. The 45 
WS, SMG and National Weather 
Service in Melbourne, Fla. (NWS 
MLB) have the responsibility of issu- 
ing weather forecasts for airfields 
located in central Florida. SMG and 
45 WS share forecasting responsibil- 
ity for the SLF depending on the mis- 
sion. The 45 WS has forecasting re- 
sponsibility for the CCAFS Skid Strip 


and Patrick Air Force Base (PAFB) 
while the NWS MLB is responsible 
for issuing terminal aerodrome fore- 
casts (TAF) for airports throughout 
central Florida. In the previous phase 
(Bauman 2009), Dr. Bauman calcu- 
lated lightning climatologies for the 
SLF and eight other airfields in cen- 
tral Florida based on a 19-year rec- 
ord of cloud-to-ground (CG) lightning 
data from the National Lightning De- 
tection Network (NLDN) for the warm 
season months of May through Sep- 
tember (1989-2007). The climatolo- 
gies included the probability of light- 
ning at 5-, 10-, 20- and 30-NM dis- 
tances from the center point of the 
runway at each site. The climatolo- 
gies were stratified by flow regimes 
with probabilities depicted at 1-, 3-, 
and 6-hour intervals. This phase up- 
dates the previous work by adding 


14 sites to the 9-site database in- 
cluding the CCAFS Skid Strip, PAFB 
and 12 commercial airports. It also 
adds three years of NLDN data re- 
sulting in a POR for the warm sea- 
son months from 1989-2010. In addi- 
tion to the flow regime stratification, 
moisture and stability stratifications 
will be added to separate more ac- 
tive from less active lighting days 
within the same flow regime. 

New Sites 

The NWS MLB requested an ad- 
ditional nine sites be added to the 
climatology to support adjacent NWS 
Weather Forecast Offices in Florida. 
Dr. Bauman agreed to add the sites 
if time permitted after he finished 
processing the first 23 sites. Dr. Bau- 
man requested the NLDN data for 
May-September 1989-2010 from Mr. 


Roeder of the 45 WS for the nine 
sites. The 14 WS prepared the 
NLDN data files and Dr. Bauman 
downloaded them from their servers. 

PWAT Stratification 

Dr. Bauman modified existing S- 
PLUS scripts to include the precipita- 
ble water (PWAT) stratification for 
each site using the PWAT thresholds 
from the closest sounding location: 
Jacksonville (JAX), Tampa (TBW), 
Miami (MFL) or CCAFS (AMU Quar- 
terly Report Q1 FY1 1). The PWAT 
stratification threshold values varied 
by up to 13% among the four sound- 
ing locations in any given warm sea- 
son month. He then generated new 

Vandenberg Air Force 
Base North Base Wind 
Study (Mr. Wheeler) 

The 30 WS states that terrain in- 
fluences along the extreme northern 
fringes of Vandenberg Air Force 
Base (VAFB) make it difficult for fore- 
casters to issue timely and accurate 
high wind warnings for that part of 
the base during northeasterly wind 
events. These events tend to occur 
during the winter or early spring 
when they are under the influence of 
the Great Basin high pressure weath- 
er regime. The LWOs have seen the- 
se rapid wind increases in Towers 
60, 70 and 71 (Figure 3) along the 
northern edge of VAFB in excess of 
the 35 kt warning threshold. For this 
task, the 30 WS requested the AMU 
analyze data from days when these 
towers reported winds in excess of 
35 kt and determine if there are any 
precursors in the observations that 
would allow the LWOs to better fore- 
cast and warn their operational cus- 
tomers of these wind events. 

VAFB Wind Tower Database 

Figure 5 is a Google Earth map 
showing the locations of all wind tow- 
ers on VAFB. Towers 60, 70 and 71 
along the northern part of VAFB are 
the primary wind towers Mr. Wheeler 
used for this study. The POR is Janu- 
ary 2004 through March 2010. During 
this period the tower data came from 
two sensor types, each formatted 


lightning climatologies for the 23 
original sites using the PWAT stratifi- 
cations. 

Because the task was ahead of 
schedule, he processed the NLDN 
data for the additional nine sites re- 
quested by NWS MLB using existing 
S-PLUS scripts. He also completed 
the climatologies for 5 sites, resulting 
in completed climatologies for 28 of 
the 32 sites. 

Upcoming Work 

Dr. Bauman will finish the light- 
ning climatologies for the remaining 
four sites and will update the graph- 
ical user interface (GUI) with the 

differently. The data from January 
2004 to October 2007 were from the 
mechanical (cup and vane) wind sen- 
sors and were archived every 5 
minutes. The data from November 
2007 through March 2010 were from 
the ultrasonic sensors and were ar- 
chived every 1 minute. 

Ms. Crawford extracted 5-minute 
peak wind speeds from the 1 -minute 
ultrasonic data using S-PLUS. This 



Figure 5. Google Earth map of the 
VAFB tower locations as yellow 
circles with white tower numbers. 
The towers for this task are 
surrounded by the yellow ellipse. 


PWAT -stratified values and addition- 
al sites. He will deliver the updated 
GUI to the customers prior to the be- 
ginning of the 2011 warm season. He 
will then assess sounding stability 
stratifications that have a relatively 
high correlation to lightning occur- 
rence as found in previous AMU 
work. This will only include stability 
parameters available to the forecast- 
ers in their weather analysis and dis- 
play systems. 

For more information contact Dr. 
Bauman at 321-853-8202 or 

bauman.bill@ensco.com. 


process was needed so all the tower 
data would have the same 5-minute 
time resolution. Mr. Wheeler import- 
ed the new 5-minute data into the 
existing 2004-2007 Excel database 
so he could complete the analysis of 
each tower’s peak wind direction and 
speed for each event. He created 
Excel charts that showed the wind 
direction and speed patterns to help 
him find precursors to high wind 
speed events. He also retrieved sur- 
face weather maps, 850 and 500 mb 
upper air maps, and VAFB sounding 
data for each of the event days from 
Plymouth State University (htl 
vortex.plvmouth.edu/) and the Na- 
tional Centers for Environmental Pre- 
diction (NCEP) ( http://www.hpc. 

i.qov/dailvwxmap/) in order 
to categorize and detail the weather 
on each event day. 

Wind Tower Events 

The 30 WS identified 66 event 
days from their cool seasons 
(October-March) in the years 2004- 
2010. Of the 66 event days, several 
did not have any tower data and oth- 
ers did not have any observations 
that met the 35 kt criterion. There 
were only 30 event days in which the 
winds at the study towers met or ex- 
ceeded 35 kt. Mr. Wheeler created 
Excel charts of wind direction and 
speed for each event day and wind 
roses for each event period. He ana- 
lyzed the surface plots to determine 
the synoptic pattern and placement 



of surface high pressure systems, and 850 and 500 
mb charts to determine the strength of synoptic sys- 
tems. He also created Google Earth maps to look for 
orographic changes that would cause a 35-kt-plus 
wind event at any of the three towers. 

Data Analysis 

Mr. Wheeler created Excel plots of wind direction 
and speed for each of the 30 event days. He noted a 
trend on most event days in which the prevailing winds 
at 54 ft would switch to the northeast prior to peak 
winds exceeding 35 kt. The timing between the switch 
and the 35-kt observations ranged from 10 minutes to 
2 hours. The key similarity was that all wind gusts 
meeting or exceeding 35 kt were always from the 
northeast on each event day. Charts of peak wind di- 
rection and speed for each event day clearly show 
that the measured peak wind speeds were from a 
northeast trajectory. 

Example Event Day 

The 30 WS identified 28 November 2007 as one of 
their strong wind event days. Synoptic conditions on 
this day (Figure 6) showed high pressure dominating 
the Pacific Northwest with a ridge extending southeast 
into Nevada. With the ridge in this location, the sur- 
face winds should have been out of the east- 
southeast at VAFB. 
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Figure 6. Surface weather map on 28 November 
2007 at 1200 UTC High Pressure in the Pacific 
Northwest with a ridge toward Nevada. The red 
star in the lower-left marks VAFB. (h[ 

www. hoc, ncep. noaa, Qov/dailvwxmap/) . 


Hourly Tower 5-min 54 ft Peak Winds 



Hourly Tower 5-min 54 ft Peak Winds 
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Figure 7. The peak wind direction (top) and speed 
(bottom) on 28 November 2007 0230-2000 UTC. The 
orange line marks the first 35 kt peak speed occurrence 
on Tower 70. 


The wind direction and speeds for this event at Tower 
70 are shown in Figure 7. This was the first northern tow- 
er to report the wind shift. The wind pattern at the three 
northern towers began with a southeast to south wind of 
around 10 kt. Around 0400 UTC; the 54-ft winds began 
backing to the northeast and increasing in speed. By 
0430 UTC the winds at Tower 70 showed a distinct north- 
east shift with an increasing wind speed. At 0500 UTC, 
winds at Tower 70 become steady out of the northeast. 

By 0530 UTC, the peak wind speed at Tower 70 was 35 
kt from the northeast. Over the next several hours the 
wind continued to gust above 40 kt. The wind shift at 
0500 UTC occurred 30 minutes before the first recorded 
35-kt peak wind. The winds weakened by 1800 UTC and 
shifted to the northwest. 

Another way to show the wind flow and magnitude of 
peak wind speeds is to use wind rose charts. Figure 8 
shows the wind rose for Tower 70 during the same time 
period as Figure 7. The red dots are wind speeds >35 kt 
at 54 ft. This chart clearly shows that all of the wind 
speeds >35 kt were from the northeast on this event day. 

During the winter months when cold high pressure 
systems build in toward northern Nevada, as shown in 
Figure 6, surface winds at VAFB would generally be out 
of the northeast to east. A Google Earth terrain map of 




period as Figure 7. The red dots represent wind speeds > 
35 kt. Other wind speed ranges and associated colors are 
in the top right legend. 



Figure 9. Google Earth terrain map showing the 
three northern VAFB tower locations as yellow 
circles with white tower numbers. The white dash 
line highlights top of range to the east and the white 
arrows highlight valleys that would enhance wind 
flow out of the northeast. 


northern VAFB (Figure 9) highlights 
a north-south ridge to the east. The 
slope of the terrain shows valleys 
oriented southwest to northeast near 
each tower leading toward the top of 
the ridge. These valleys likely help 
focus and strengthen a northeast to 
east surface wind. 

Recommendations 

A northeast tower wind flow pat- 
tern was apparent on all of the 30 
event days. Mr. Wheeler showed 


through his analysis that the cases in 
this task all showed an apparent shift 
of the surface winds to a northeast 
flow pattern 20 minutes to 2 hours 
before the winds reached 35 kt. The 
next step in developing a high-wind 
alert capability for VAFB would be 
the installation of a local mesoscale 
model that would incorporate all of 
the local data sets. The model 
should then be tested to see if it can 
forecast these types of events with a 
lead time of 2 to 24 hours. This 


would allow the meteorologist at 
VAFB to alert the operational cus- 
tomers in a timelier manner so pro- 
tective action could be taken. 

Status 

Mr. Wheeler completed his re- 
view and analysis of the VAFB wind 
events and started writing the final 
report. 

For more information contact Mr. 
Wheeler at 321-853-8205 or 

wheeler.mark@ensco.com 
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INSTRUMENTATION AND MEASUREMENT 


Upgrade Summer 
Severe Weather Tool 
Phase III (Dr. Watson) 

The 45 WS Commander’s morn- 
ing weather briefing includes an as- 
sessment of the likelihood of local 
convective severe weather for the 
day. This forecast is provided in or- 
der to enhance protection of person- 
nel and material assets of the 45th 
Space Wing, CCAFS, and KSC. The 
severe weather elements produced 
by thunderstorms include tornadoes, 
convective surface winds of 50 knots, 
and/or hail with a diameter of 0.75 
inches. Forecasting the occurrence 
and timing of these phenomena dur- 
ing the warm season (May- 
September) is challenging for 45 WS 
operational personnel. In previous 
tasks, the AMU analyzed stability pa- 
rameters and synoptic patterns from 
Central Florida severe weather days 
during 1989-2003 to determine which 
were important to severe weather 
development (Bauman et al. 2005). 
The AMU then created a hyper text 
markup language (HTML) tool using 
the important parameters and pat- 
terns to help determine the probabil- 
ity of issuing severe weather watches 
and warnings for the day. The HTML 
tool was replaced with a Meteorologi- 
cal Interactive Data Display System 
(MIDDS) GUI in a follow-on task 
(Wheeler 2009) that retrieved stability 
parameters and other information 


from MIDDS automatically, minimiz- 
ing the forecaster’s interaction with 
the tool. Later, the AMU updated the 
severe weather database with data 
from the years 2004-2009, re- 
analyzed the data to determine the 
important parameters, made appro- 
priate adjustments to the index 
weights depending on the results of 
the analysis, and updated the MIDDS 
GUI (Wheeler 2010). For this task, 
the 45 WS requested the AMU up- 
grade the severe weather database 
by adding weather observations from 
2010, update the verification data set 
with results from the summer of 
2010, use statistical logistic regres- 
sion analysis on the database and 
develop a new forecast tool if appro- 
priate, and update the MIDDS GUI, if 
necessary. 

Logistic Regression Analysis 

Following Lambert and Wheeler 
(2005), Dr. Watson conducted pre- 
dictor selection using S-PLUS, which 
has a built-in logistic regression func- 
tion. She used a process called 
screening regression to determine 
which candidate predictors to include 
in the logistic regression equation. In 
this approach, predictors were added 
to the equation one at a time. At each 
step, the candidate predictor that cre- 
ated the biggest reduction in the re- 
sidual deviance was chosen as the 
next predictor in the equation. Figure 
10 shows the percent reduction in 
residual deviance from the NULL 

model as each 
predictor was 
added. The To- 
tal Totals (TT) 
reduced the 
residual devi- 
ance by the 
most (~8%) 
and was cho- 
sen as the first 
predictor in the 
equation. The 
second predic- 
tor was the flow 
regime, which 
brought the to- 
tal reduction of 


residual deviance to ~13%. The Lift- 
ed Index (LI) and upper-level jet ex- 
istence were the third and fourth pre- 
dictors in the equation, respectively, 
producing the final reduction in resid- 
ual deviance of 15%. 

Logistic Regression Equation 
Performance 

Dr. Watson compared the calcu- 
lated forecast probabilities using the 
four predictors in Figure 10 and de- 
termined their performance using the 
binary severe weather observations 
and four tests that measure forecast 
performance. The tests included 

• Mean Squared Error (MSE), 

• Brier Skill Score (BSS), 

• Distributions of the probability 

forecasts for days with and with- 
out severe weather, and 

• Contingency table statistics. 

The MSE is the mean of the 
squared differences between the 
forecast probabilities and the obser- 
vations. The MSE for a perfect fore- 
cast is 0, with larger MSE indicating 
decreasing accuracy of the forecast. 
The MSE was computed for the lo- 
gistic regression equation using the 
development and verification da- 
tasets. The MSE for the full develop- 
ment dataset was 0.10, which indi- 
cates skill in predicting severe weath- 
er. However, when the data were 
split into severe and non-severe 
events, the MSE was 0.61 and 0.03, 
respectively. Similarly, the MSE for 
the full verification dataset was 0.1 1 , 
but 0.59 and 0.02 for severe and non 
-severe events, respectively. These 
results indicate that the equation was 
biased towards predicting non-events 
and failed to adequately predict se- 
vere weather events. The MSE was 
also computed for the Total Threat 
Score (TTS) using the verification 
dataset. The MSE for the full verifica- 
tion dataset was 0.07 and was 0.26 
and 0.03 for severe and non-severe 
events, respectively. Based solely on 
MSE, the TTS was a better predictor 
of severe weather events than the 
logistic regression equation. 



Figure 1 0. The total percent reduction in residual 
deviance from the NULL model as each predictor was 
added to the equation using the development dataset. 


Forecast Probability Distributions 
for Severe and Non-Severe Weather Days 



Table 2. Skill scores for the TTS and logistic re- 
gression forecasts. 

Skill Score 

TTS 

Logistic Regression 

POD 

0.73 

0.35 

FAR 

0.23 

0.42 

CSI 

0.60 

0.28 

HSS 

0.70 

0.36 

TSS 

0.68 

0.30 


Figure 11. Forecast probability distributions for severe 
(red) and non-severe (blue) days in the verification data. 
The y-axis values are the frequency of occurrence of 
each probability value, and the x-axis values are the 
forecast probability values output by the equation. 


The BSS measures the differ- 
ence in skill of the logistic regression 
equation against a reference fore- 
cast. The BSS denotes a percent im- 
provement (degradation) in skill of 
the equation over the reference fore- 
cast when it is positive (negative). 
The TTS for the verification dataset 
was used as the reference forecast. 
The BSS values for the verification 
dataset were -57% for the full da- 
taset, -131% for the severe weather 
events, and 34% for non-severe 
weather. As with the MSE, these re- 
sults indicate that the logistic regres- 
sion equation is biased towards pre- 
dicting non-events as the percent 
improvement for the non-severe 
weather is large. However, the per- 
cent degradation for predicting se- 
vere events is quite large, again indi- 
cating that TTS is a better tool for 
predicting severe weather. 

The equation probability fore- 
casts from the verification dataset 
were stratified by severe and non- 


severe weather 
days. The dis- 
tribution of the 
probability val- 
ues was calcu- 
lated for each 
stratification. 
Figure 1 1 shows the probability distri- 
bution for severe days (red curve) 
and non-severe days (blue curve). If 
the equation performed well, the red 
(blue) curve would have a minimum 
(maximum) in the lower probability 
values that increase to a maximum 
(minimum) at the higher values. The 
non-severe weather days had a peak 
frequency near 65% at probability 
values of 0.1 and then decreased to 
near 0 at 0.6. It shows a high per- 
centage of low probabilities for non- 
severe events and a low percentage 
of high probabilities as expected for 
good performance. The severe 
weather days had a small peak of 
30% at probability values near 0.2 
followed by a dip and then another 
small peak near 15% at probability 
values at 0.4. This indicates that the 
equation performed poorly for severe 
weather days. The maximum at 0.2 
and minimum at 0.6 suggests the 
equation is under-forecasting severe 
weather events. 


Table 2 shows the contingency 
table statistics for the TTS and lo- 
gistic regression equation probabili- 
ties for the verification dataset. The 
Probability of Detection (POD) and 
Critical Success Index (CSI) are 1 for 
a perfect forecast and 0 for no skill, 
and vice versa for the False Alarm 
Rate (FAR). The Heidke Skill Score 
(HSS) and True Skill Statistic (TSS) 
are 1 for a perfect forecast, 0 for per- 
formance equal to a random forecast, 
and < 0 for performance worse than 
that of a random forecast. It is evi- 
dent that the TTS outperforms the 
equation in every computed statistic. 

Final Report 

Dr. Watson completed the first 
draft of the final report. She then 
modified the report based on recom- 
mendations received from the inter- 
nal AMU and external customer re- 
views. She submitted it for NASA ap- 
proval and will have the report up- 
loaded to the AMU website after that 
approval is received. 

For more information contact Dr. 
Watson at 321-853-8264 or 

watson.leela@ensco.com 
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MESOSCALE MODELING 


MesoNAM Verification 
Phase II (Dr. Watson) 

The 45 WS LWOs use the 12-km 
resolution North American Mesoscale 
(NAM) model (MesoNAM) text and 
graphical product forecasts exten- 
sively to support launch weather op- 
erations. In Phase I of this task 
(Bauman 2010), the AMU measured 
the actual performance of the model 
objectively by conducting a detailed 
statistical analysis of model output 
compared to observed values. The 
model products included hourly fore- 
casts from 0 to 84 hours based on 
model initialization times of 00, 06. 

12 and 18 UTC. The objective analy- 
sis compared 3.5 years of MesoNAM 
forecast winds, temperature and 
dewpoint, as well as the changes in 
these parameters over time, to the 
observed values from the sensors in 
the KSC/CCAFS wind tower network. 


For this task, the 45 WS requested 
the AMU modify the current tool by 
adding an additional year of model 
output to the database and recalcu- 
lating the verification statistics. The 
AMU will also update the current GUI 
with the new statistics. This tool 
helps the LWOs understand the mod- 
el’s performance when they use it to 
evaluate LCC during launch opera- 
tions. 

Wind Tower Data 

Dr. Watson acquired the KSC / 
CCAFS wind tower data for the peri- 
od February 201 0 to February 201 1 
from the AMU archive, and used the 
AMU wind tower QC software to re- 
move erroneous observations from 
the dataset. She used S-PLUS 
scripts written by Ms. Crawford to 
import and modify the QC’d wind 
tower observation files to remove un- 
needed time periods from the dataset 


for each tower and to fill in missing 
values with the appropriate designa- 
tion. The locations of the towers used 
for the verification are shown on the 
map of KSC/CCAFS in Figure 12. 
Next. Dr. Watson modified a previ- 
ously written S-PLUS script with the 
help of Ms. Crawford to compute the 
mean value for each observed pa- 
rameter at the top of every hour us- 
ing the observations from 30 minutes 
prior and 25 minutes after the hour. 

MesoNAM Forecast Products 

Dr. Watson requested and ob- 
tained the NCEP MesoNAM fore- 
casts from Mr. Randy Nyman of AC- 
TA, Inc. She imported the forecast 
files using S-PLUS scripts. 

For more information contact Dr. 
Watson at 321-853-8264 or 

watson.leela@ensco.com 



Figure 12. Map of KSC/CCAFS showing the wind tower locations as red pentagons labeled with tower 
number and the supported launch activity, MesoNAM model grid points as green circles, and the 
CCAFS weather station as a magenta square labeled KXMR (Bauman 2010, Figure 1). 
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AMU ACTIVITIES 


AMU Chiefs Technical 
Activities (Dr. Merceret, 
Dr. Huddleston) 

Dr. Merceret began investigating 
errors in radar measurements of 
cloud location, height and thickness 
after concerns were raised by the 
Lightning Advisory Panel. He is ex- 
amining error contributions from ef- 
fects of propagation and beam ge- 
ometry as well as the computational 
processes used by the radar product 
generation software. He requested 
climatological data on microwave 
index of refraction gradients near 
KSC/CCAFS and software to exam- 
ine the propagation effects of those 
gradients. Dr. Merceret also began a 


literature search on the effects of 
beam geometry. 

Dr. Merceret’s work focused on 
effects of microwave propagation 
conditions on height measurements. 
This issue is also a concern for radi- 
ological contingency planning for the 
Mars Science Lander (MSL) sched- 
uled for launch in November 201 1 . 
He generated a variety of soundings 
of the refractivity, N. based on clima- 
tological statistics of N at KSC / 
CCAFS provided by the Air Force. 

He analyzed these soundings using 
the ray tracing capability in a soft- 
ware package called Advanced Re- 
fractive Effects Prediction System 
(AREPS) obtained from the US Na- 
vy. Dr. Merceret’s preliminary exami- 
nation suggests that refractive 


height errors are unlikely to be of 
concern for evaluation of the light- 
ning launch commit criteria, but 
probably will be of concern for MSL 
contingency operations. 

In March, Dr. Huddleston took 
over the position of AMU chief from 
Dr. Merceret. Dr. Huddleston contin- 
ued to familiarize herself with the 
needs of the 45 WS and the activi- 
ties of the AMU. She corrected the 
45 WS lightning strike location anal- 
ysis software to account for an error 
that occurred if a lightning stroke 
occurred at precisely the same lati- 
tude and longitude as the center of 
the area of interest. Dr. Huddleston 
is working on a draft paper of the 
lightning probability algorithm for the 
Journal of Spacecraft and Rockets. 
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14 WS 14th Weather Squadron 

30 SW 30th Space Wing 

30 WS 30th Weather Squadron 

45 RMS 45th Range Management Squadron 

45 OG 45th Operations Group 

45 SW 45th Space Wing 

45 SW/SE 45th Space Wing/Range Safety 

45 WS 45th Weather Squadron 

AFSPC Air Force Space Command 

AFWA Air Force Weather Agency 

AMU Applied Meteorology Unit 

BSS Brier Skill Score 

CCAFS Cape Canaveral Air Force Station 

CG Cloud-to-Ground lightning 

CSI Critical Success Index 

CSR Computer Sciences Raytheon 

FAR False Alarm Rate 

FSU Florida State University 

FY Fiscal Year 

GSD Global Systems Division 

GUI Graphical User Interface 

HSS Heidke Skill Score 

HTML Hyper Text Markup Language 

JAX Jacksonville, Fla. 3-letter identifier 

JSC Johnson Space Center 

KSC Kennedy Space Center 

LCC Launch Commit Criteria 

LI Lifted Index 


LWO 

Launch Weather Officer 

MesoNAM 12-km North American Mesoscale model 

MFL 

Miami, Fla. 3-letter identifier 

MIDDS 

Meteorological Interactive Data Display 
System 

MSE 

Mean Square Error 

MSFC 

Marshall Space Flight Center 

NCEP 

National Centers for Environmental 
Prediction 

NLDN 

National Lightning Detection Network 

NOAA 

National Oceanic and Atmospheric 
Administration 

NWS MLB National Weather Service in Melbourne, 1 

PAFB 

Patrick Air Force Base 

POD 

Probability of Detection 

POR 

Period of Record 

PWAT 

Precipitable Water 

QC 

Quality Control 

SLF 

Shuttle Landing Facility 

SMC 

Space and Missile Center 

SMG 

Spaceflight Meteorology Group 

TAF 

Terminal Aerodrome Forecast 

TBW 

Tampa, Fla. 3-letter identifier 

TSS 

True Skill Statistic 

TT 

Total Totals 

TTS 

Total Threat Score 

USAF 

United States Air Force 

VAFB 

Vandenberg Air Force Base 
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The AMU has been in operation since September 1991. Tasking is 
determined annually with reviews at least semi-annually. 

AMU Quarterly Reports are available on the Internet at http://science.ksc.nasa.gov/amu/ . 

They are also available in electronic format via email. If you would like to be added to the email distribution list, 
please contact Ms. Winifred Crawford (321-853-8130, crawford. winnie@ensco. co r). 

If your mailing information changes or if you would like to be removed from the distribution list, please notify 
Ms. Crawford or Dr. Lisa Huddleston (321-861-4952, Lisa.L.Huddleston@nasa.gov) . 
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